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Abstract 

With this paper we begin an investigation of difference schemes 
that possess Darboux transformations and can be regarded as natural 
discretizations of eUiptic partial differential equations. We construct, 
in particular, the Darboux transformations for the general self ad- 
joint schemes with five and seven neighbouring points. We also intro- 
duce a distinguished discretization of the two-dimensional stationary 
Schrodinger equation, described by a 5-point difference scheme involv- 
ing two potentials, which admits a Darboux transformation. 
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1 Introduction 



Linear differential operators admitting Darboux Transformations (DTs) play 
a crucial role in the theory of integrable systems IH I21- i) One can asso- 
ciate with such operators integrable nonlinear partial differential equations 
(PDEs). ii) One can make use of the spectral theory of such linear oper- 
ators to solve classical initial - boundary value problem for the nonlinear 
PDEs. iii) One can construct solutions of the nonlinear equations from sim- 
pler solutions, using the Darboux - Backlund transformations, iv) One can 
often associate with the nonlinear PDE a geometric meaning, inherited by 
the geometric properties of the linear operator. 

It is natural to search for a discretization of this beautiful picture. In 
general, searching for distinguished discretizations of linear differential oper- 
ators that admit general DTs (integrable discretizations) is not a trivial task. 
Several methods of "integrable discretization" have been used so far (see e.g. 
jni), but no one is fully satisfactory. Surprisingly enough, a distinguished inte- 
grable discretization of the two-dimensional stationary Schrodinger equation 
is, to the best of our knowledge, not present in the literature and one of our 
goals is to change this situation. 

In recent years the study of linear difference equations that admit DTs 
were undertaken. Most of the results are based on the 4-point difference 
scheme, i.e. a scheme that relates four neighbouring points 

m,n+l ~l~ ^m,n'4^m,n (l) 

(where a, /3, 7 are real functions of two discrete variables {m,n) e 1? and 
where the standard notation fm,n = f{^,n) is used throughout the paper), 
which is proper for discretizing second order hyperbolic equations in the 
canonical form. In particular, the proper discrete analogue of the Laplace 
equation for conjugate nets 

^,u,+C^,u+D^,v=0, (2) 

whose DTs were obtained in ^U], turns out to be the 4-point scheme jH] 

'^m+l,n+l Olm,n'4^m+l,n ~l~ Pm,n'4^m,n+1 ~l~ (1 Oim,n Pm,n)'^m,ny (3) 

describing a lattice with planar quadrilaterals [7j, [S], whose general DTs 

were extensively studied in recent years (see, e.g., [HI CHI El)- While the 
Moutard equation [T2] 

*,™=i^^, (4) 
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relevant in the description of asymptotic nets, splits naturally in the discrete 
case into two equations jlHl HI] 




(5) 



admitting DTs. Laplace transformations for the general 4-point scheme 
are also known [71 UHl 11^] • 

The discretization of elliptic equations in two dimensions admitting gen- 
eral DTs is much less studied and understood. In the continuous case, elliptic 
equations and their Darboux transformations are often obtained from the hy- 
perbolic ones regarding {u, v) as complex variables: 



but this approach presents some problems in the discrete case. 

In our opinion, a proper discretization of an elliptic operator should satisfy 
two basic properties. 

• It should be applicable to solve generic Dirichlet boundary value prob- 
lems on a 2D lattice; 

• It should possess a class of DTs which must be (at least) as rich as that 
of its differential counterpart. 

It is very easy to convince one-self that the first criterion cannot be sat- 
isfied by the 4-point scheme (P) and that at least 5-point difference schemes 
should be introduced [T7j. The most general 5-point scheme reads as follows: 

0'm,n'4^m+l,n ~l~ ('^m— l,n ~l~ '^m,n)'^m—l,n ~l~ ^m,n'^m,n+l ~l~ (^m,n— 1 ~l~ ^m,n)'^m,n—l 



where a, b, w, z are functions of the discrete variables (m,n) G Z^. 



u = X + iy 



V = X — ly, 



(6) 




(7) 
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(m,n+l) 



(m-l,n) 



(m,n) 



(m+l,n) 



' (m,n-l) 

Fig.l The 5 - point stencil 
Its natural continuous limit 

(em, en) {x,y), (8) 
in which the lattice spacing e goes to zero and 

hm,n ) \ B{x,y) I 

(9) 

fm,n = 2{A + B)-t{A,., +B,y +W + Z)+e^[F{x, + +B,yy )] +0(e3), 

gives the following second order PDE in the independent variables (x, y) 

{A^,, +W^„ ),y +Z^<,y = F^, (10) 

in which the ^ ,xy term is missing (we remark that the ,xy term is instead 
the only second order term present in equations and (jH)). Equation (fTUj) 
is elliptic if AB > and its canonical form is obtained by setting A = 1 = B. 
The functions w, z, W, Z measure the departure of equations ((Tj) and (fTUIl 
from self-adjointness, which is a basic property in most of the applications. 
In this paper we present the following results. 

• We construct the DTs for the self-adjoint reductions of the 5-point 
scheme d?)) and of its continuous limit (fTUj) . 

• We use the gauge covariance property of the general self-adjoint 5-point 
scheme to construct a distinguished integrable discrete analogue of the 
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stationary Schrodinger operator in two dimensions. "Integrable" in the 
sense that, not only it reduces to the stationary Schrodinger operator 
in the continuous hmit, but also it possesses DTs. 

• We show that the construction of DTs for the self-adjoint 5-point 
scheme can be applied to the case of self-adjoint schemes involving 
more neighbouring points, illustrating such a generalization in the case 
of the general self-adjoint 7-point scheme, for which Laplace transfor- 
mations are already known ^Sl ^| • 

The paper is organized as follows. Section |21 is devoted to the construction 
of the DTs for the most general self-adjoint 5-point difference operator with 
complex coefficients: 

where and T„ are the translation operators with respect to the discrete 
variables {m,n) 6 "I?: 

Tmfm,n fm+l,n) Tnfm,n fm,n+l 

and a, b, c are complex functions of the two discrete variables {m,n) G Z^. 
We also use the following difference operators 

^mfm,n fm+l,n fm,n ^nfm,n fm,n+l fm,n 

^—mfm,n fm—l,n fm,n ^—nfm,n fm,n—l fm,n- 

The operator £5 is formally self-adjoint with respect to the bilinear form: 

< f,9 >■= X] fm,n9m,n (12) 
m,n 

and, in the particular case a = 6 = 1, it reduces to the 5-point scheme 

Csch = Tm + + T„ + - /, (13) 

which is the simplest and most used (in numerical applications) discretization 
of the Schrodinger operator 

+ (14) 



while the further constraint / = 4 leads to 



C^ = T^ + T^^ + r„ + T-^ - 4, (15) 

i.e. to a discretization of the Laplacian ^ + ^ often used in the numerical 
studies of elliptic boundary value problems [EI . 

In the continuous limit, the operator £5 goes to the second order partial 
differential operator 

L = Ai, + B^ + A,^j, + B,yl-F (16) 

which is formally self-adjoint with respect to the bilinear form: 

< fi9 >■= / fgdxdy. (17) 



We notice that, while both operators (|TT|l . ([THjl are self-adjoint, the two 
"integrable" discretizations (0) of the self-adjoint Moutard equation (jH) are 
just mutually adjoint. 

In section IHl using the following covariance property (gauge invariance) 
of the operator £5: 

^m,n ^ ^m,n ^m,n9m,n9m+l,n^ ^m,?i ^ ^m.n ^m,,n9in,n9in,n+li (-^^) 

fm,n ^ fm,n fm,n9r, 



we draw the operator £5 to the form 

/I ^m,n rp I rm-ljnrp-l , ^m,n rp , Tmn-i -j^ /1n^ 

-t-ScWni — + 1^ +- ln + —^ J„ - /• \^^) 

J- m+l,n ^ m,n ^ m,n+l ^ m,n 

The operator (fTn|) . symmetric with respect to 7r/2 rotations in the plane grid 
{■m,n), reduces to the two-dimensional Schrodinger operator (fT^ under the 
natural continuous limit and possesses DTs. Therefore it appears as the "in- 
tegrable" discrete analogue of the two-dimensional Schrodinger operator. It 
is therefore the proper starting point in the search for the integrable discrete 
analogues of the nonlinear PDEs of elliptic type associated with ()14j) (like 
the Veselov-Novikov hierarchy nonlinear a models 03 ^^^1 the Ernst 
equation 
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The corresponding gauge transformation for the operator L in is 

L->L = gLg, 
A = g'A, B = g'B, (20) 
F = gLg. 

We draw the reader's attention to the fact that constraint A = B is 
conserved by the gauge transformation ()20j) and the same is true for the 
constraint am,n+iflm,n = bm+i,nbm,n with respect to gauge (fTHj). In the case 
of the reduction A = B^ the operator ()16|) can be gauged into the form of 
the 2D Schrodinger operator (fT^ while, in the case am,n+iOm,n = &m+i,n&m,n, 
the operator (fTTj) can be gauged into the operator (fT5|) . But while there 
exist DTs which preserves the form of equation (fT^ . we don't know DTs 
that preserves the form of equation (fT^ . This is the reason why we view 
the operator (fT^ and not the operator (fT^ as the integrable discretization 
of the Schrodinger equation (fT^ . 

Finally, in section |^ we show that there is another self-adjoint scheme 
involving more neighbouring points which possess DTs of the type presented 
in section 121 It is a self-adjoint 7-point difference scheme 



bm,n—l'^m,n—l ~l~ Sfn+l,n'4^m+l,n—l ~l~ '^m,n+l'^m—l,n+l fm,n4^ri 



(21) 



which, in the continuous limit 



\ 



J 



A{x,y) 
B{x,y) 

V -S{x,y) ) 



+ 0(6^ 



(22) 
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= 2(A+5+^)-e(A„. +5„ +2^„. +25,, )^e\F[x, i/) + -(A„., +5,,^, )\^0[^) 

goes to the most general second order self-adjoint differential equation in two 
independent variables: 



(A^„ )„ +(5^,J, ),, +(5^,y ),y +(5^„ 



(23) 



The operator £7 in (j^U, introduced in [TH], admits the representation 
£7 = QQ~^ + w in terms of a 3-point difference operator Q and of its adjoint 
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Q"*", and this factorization plays a crucial role in the construction of Laplace 
transformations for £7 jTHl UH] and in the development of an associated dis- 
crete complex function theory [22] ■ The 5-point operator £5 does not admit 
the above representation and the construction of the DTs of this paper does 
not make use of it. We also remark that several discretizations of the ID 
Schrodinger operator have appeared in the literature throughout the years 
(see, f.i., 12211211123 and [1^] with references therein included). We end this 
introduction observing that a 4-point scheme with a complexification of the 
discrete variables (m, n) analogous to ©, was used in [2^1 to construct a 
discrete analogue of the a model. 

2 5-point self adjoint operator and its Dar- 
boux transformation 

In this section we present a Darboux transformation for the equation: 
where a^^n, hm,n and fm,n are given functions. Let 6 be another solution of 



then: 

fm,n '0 iflm,n^m+l,n ~l~ (^m—l,n^m—l,n ~l~ ^m,n^in,n+l ~l~ ^m,n—l^m,n—\) • (^6) 

Eliminating fm,n from ()24|) and ()25|) we get 

^miflm—l,n'4^m,n^m~l,n ^m—l,n^m,n'4^m—l,n)~\~ 
^n(fim,n—l'4^m,n^m,n—l ^m,ri~l^m,'n^Pm,n—\) fl- 
it means that there exists a function a such that 



(27) 



A n/ — —h n n A 'Vrn,n \'^° ) 

'~'m,n—l"m,n"m'n—l'-^—n a 
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Setting 

, / 0(m,n 

rm,n n 

"ra,n 

we find that „ satisfies the following equation 
where 

^m-l,n ~ 7 7^ "m,n-l ~ B V'^^/ 

"m— l,n— 1 "»n— 1 '^m— l,n— 1 "m— l,n— 1 

and 

\ "m+l,n "m— l,n "m,n+l "m,n—l / 

(31) 

Comparing equations and (jHH), we also infer that 6*' = 1/6* is a solution 

of (121. 

In the continuous limit with 6'm„ = 6(x, ?/) + O(e^) and with 
expanded according to to get (jHt"), we obtain the DT 

i(evi>')„ = -Ae(|),. ^(ev^O,. = Be{^u (32) 

from the solution space of equation 

(A^,,,),,+(5^„)„=F^ (33) 
to the solution space of equation 

(A'vl/'„.),.+(i?'^'„)„=i^'^', (34) 
where 6 is another solution of so that 

F = i[(Ae,j,,+(5e„)„], (35) 

and the new "potentials" A', B' and F' are related to the old ones as follows 

4' — J- R' — 1 
F' = 0[(A'(^),J,.+(5'(^)„)„]. ^"'^^ 

Comparing equations (jH^ and (jHSfc), we also infer that 9' = 1/0 is a solution 
of (jSl. 
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3 A two-dimensional Schrodinger operator and 
its Darboux transformation 



m,n 



The operator (fTTj) can be gauged into the form (fT^ . Indeed, if a^.n and h 
are given functions, one can find a function „ such that = ^""'"^^ and, 

under the gauge (fTH|) . we get dm,n = &m,n =: Finally, the operator C 

given by £5 = y^^C^y^^ is of the wanted form (fT^ . 

Combining the DT of the previous section with the above gauge trans- 
formation, one obtains the following DT for the discrete analogue ()19|) of the 
Schrodinger operator. 

Let Njn^n be a solution of the integrable discrete analogue of the 2D 
Schrodinger operator: 

r77i,n AT" I Tm—l.n J\T 1 Tm.n AT" 1 '^m.n — 1 AT" AT" 

r ~ -^*m+l,n ~r -/V772— l.n ~r p ~"-^*m,,nH-l ~r -^*m,ri— 1 Jm.n-^^m.n? 

-L m+l,n -L m.ri ' J- m,n+l -L m.Ti 

(37) 

and let be another solution of it. It means that 

Jm,n ~Q I '^m+l,n ' -p "m— l,n ~r z:; "r7i,n+l ~r _ f^m,7i— 1 

" VJ- m+l,n J- m,n J- m,n+l J- r7i,n 

(38) 

Eliminating /m,n from (jHHj) and (j?fj) we get that there exists a function A^^ m 
given in quadratures by 



'^m,n"m,n— 1 I To To 7 / 

J- m,n \ ^m.n ^ m.n — l J 

A /I'm 

I p J- m,n-^^ m,n I ^m,n^m—l,n \ To To ; / 



(39) 



which satisfies 



r' r' r' r' 

r' m+l,n ' r' m— 1,71 ' r' , , m,n+l r' m,n—l J m,n m,ni 

(40) 

where 



r'2 = 








f = 
J m,n 


r 


r' 

m,n — l ( 


r 


v ^ 

m.n 





Jm,n~ r I r' , , lr'i?''™+i'" r' lr'i?''"^~i'" r' , , vr'i?''"^''^+i~'~ 



^ )m,n—l\ 

(41) 
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Now the function {}' = is a solution of (jin|) . In the continuous hmit, with 
Tm.n = G{x,y) + 0(e) and with fm,n expanded according to equation (jSHl), 
one obtains the classical Moutard transformation 

= e (I) ,y , ^{QN'),y = -e (I) (42) 

from the solution space of equation 

= (43) 

to the solution space of equation 

where B is another solution of PHj) . so that 



1 

e 

and 



F ^ (0)XX ~I~0)3/J/ ) ; 



1 

ey — ye 

Furthermore 0' = 1/0 is a solution of 



4 7-point self-adjoint operator and its Dar- 
boux transformations 

We end this paper showing that the construction of DTs presented in the 
previous two sections applies also to a self-adjoint scheme involving more 
than 5 points, namely in the case of the self-adjoint 7-point scheme: 

(^m,n'^m+l,n ~l~ (^m—l,n'^m—l,n ~l~ ^m/nV^m/ri+l ~l~ ^m,n— (45) 
"5m+l,n V'm+l,n— 1 ~l~ Sm,n+l'^m—l,n+l fm,n'4^in,n) 

where a^.n, ^m.n, Sm,n and fm,n are given functions, which, as we have seen 
in the introduction, is a discretization of the most general second order, self- 
adjoint, linear, differential equation in two independent variables. 
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Let 9m,n be another solution of equation (^3]) : 

(^m,n(^m+l,n '^m— l,n^m— l,n ~l~ ^m,,n^m,n+l ~l~ ^m.n— l^m.n— (46) 
^m+l,n9m+l,n—l '5m,n+l^m— l,n+l fm,n9m,n- 

Eliminating fm,n from (ji^j) and pUj) we get 

m[^m—l,n"m,n"m—l,n\n a , j ~r "Jm/raCm— l,n"m,i 



/'4>n 





'i/'m— l.n 




^m — l,n 


V^m — 1 ,ri 


V'm,n-1 > 


— 1 , 


Om.n—1 ' 



(47) 



It means that there exists a function a such that 



^mC^m,n (^m,n— l^m,n^m,n— 1 ~l~ ^m,rfim—l,rfim,n—l)^—n~a l~ '^m,n9m—l,n9m,n— 

(48) 

Introducing 

we find that ip'm n satisfies the following equation 

a' ol/ _|_ f,l r,!/ _|_ y r.!/ I IJ r.!/ _|_ 

m,n r m+l,n m— l,n rm— l,n ' m,n r m,n+l ' m,n—l r m,n— 1 ' rzlQ"! 

s' ?// + ?// = f ih' ^ ^ 

■^m+ljn r m+l,n— 1 ' '^m,n+lTm—l,n+l Jm,nTm 



where the new fields are given by 

/ , n 4- 1 , n '^m — 1 , n 



a, 



"^i" 6m,n — lPm,n 



,/ Sm — l,n—l^m—l,nf^Tn.n — l 

'?m-l,n-iPm-l,n-l ' 
/ 1 



f = ft ( n' 1 _|_ ^1 1 _|_ IJ 1 _|_ JJ 1 I 

(50) 

and where Pm,n 9m,n(^m—l,nbm,n—l ~l~ ^m— l,n'5m,n'^m— l,n ~l~ ^m,n9m,n—lbm,n—l- 

Again 6'^^ = l/6m,n is a solution of (jl9|l . 

In the continuous limit (j221), with 6m,n = 0(a^,Z/) + ^(e^) and with fm,n 
expanded according to to obtain (j22I), "we obtain the DT 

^(0vi/')„=-Ae(|),.-5e(|)„ 

i(0vi/'),. = i?e(|)„+^0(|),. ^"'^ 
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from the solution space of equation 

(A^,, )„ +iS^,y +(^^,. = (52) 

to the solution space of equation 

where 9 is fixed solution of (jK^ : so 

F = ^[(Ae,, ),, +{Be,y ),y +{SQ,y +{SQ,, ),y ] (54) 

and the new "potentials" A', B', S' and F' are related to the old ones as 
follows 

4/ A f)/ B q/ 5 



AB-S'^ ' AB-S^ ' A_B-52 ' f'ccN 

F' = e[(A'(i),. ),. +{B'{^),y ),y +{S'{^),.),y + {S 'i^) ,y ),.]■ ^ 

Again 9' = 1/6 is a solution of (fSHj). 
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